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Abstract. We discuss generation of non-Gaussianity in density perturbation through 
the super-horizon evolution during inflation by using the so-called 5N formalism. 
We first provide a general formula for the non-linearity parameter generated during 
inflation. We find that it is proportional to the slow-roll parameters, multiplied by 

the model dependent factors that may enhance the non-Gaussianity to the observable 
ranges. Then we discuss three typical examples to illustrate how difficult to generate 
sizable non-Gaussianity through the super-horizon evolution during inflation. First 
example is the double inflation model, which shows that temporal violation of slow 
roll conditions is not enough for the generation of non-Gaussianity. Second example 
is the ordinary hybrid inflation model, which illustrates the importance of taking into 
account perturbations on small scales. Finally, we discuss Kadota-Stewart model. This 
model gives an example in which we have to choose rather unnatural initial conditions 
even if large non-Gaussianity can be generated. 
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1. Introduction: Success of inflation 



Inflation solved the various cosmological problems like horizon, homogeneity, flatness 
problem, just by assuming that the potential energy of scalar fields dominates the 
expansion of the universe. Moreover, inflation scenario naturally explains the origin 
of the almost scale invariant density perturbation. The observations in the past decades 
verified many of the predictions of inflation. Almost scale invariant spectrum is now 
confirmed and now we even know that the violation of the exact scale invariance is 4% 
level [Tj, which is also consistent with the simplest slow roll inflation scenario. 



1.1. Density perturbation 

The amplitude of quantum fluctuation of inflaton (j) during slow-roll inflation is 
determined by the unique relevant mass scale at that time, e.g. the Hubble rate H: 

H = ^. (1) 

However, on scales as large as the Horizon radius the meaning of the amplitude of 
field perturbation becomes subtle because it depends on the choice of gauge. For the 
flat slicing gauge, in which the trace of the spatial curvature is kept unperturbed, the 
perturbation equation for a scalar field becomes very simple and looks very similar 
to the one without gravitational perturbation [2]. Therefore the amplitude mentioned 
above can be in fact understood as that in the fiat slicing gauge. This amplitude 
of perturbation can be interpreted as the dimensionless curvature perturbation on a 
uniform energy density surface (. Transformation law is given by 

C = H5t = H^-i = K (2) 
27r0 

where 6t is the shift in time coordinate for this transformation, which is given by S4>/<p 
using the time derivative of the background scalar field. Writing the perturbation 
amplitude in terms of ( is useful because the curvature perturbation does not evolve on 
super-horizon scales when the evolution path of the universe is unique. This constancy 
of ( simplifies the analysis of density perturbation during infiation for a single infiaton 
model. An extension of this argument to the case of multi-component inflaton is the 
SN formahsm [31 IH El El El E] , as we explain below. 



1.2. Further steps from observations 

Further detailed comparison between the theoretical predictions and future observations 
is awaited. One important observation will be tensor-type perturbation in cosmic 
microwave background. Tensor-type perturbation is the transverse-traceless part of 
the spatial metric perturbation, which is generated by the same mechanism as the 
perturbation of the inflaton fleld. Gravitational action has a overall prefactor m^^, and 
hence we deflne the canonically normalised metric perturbation ■j/'^j, = nipiSgfj^t, so as 
to absorb this factor from the quadratic action. Then, the quadratic action for the 
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transverse traceless part of ifj^,^ becomes identical to the one for a massless scalar field. 
Therefore the amplitude of perturbation ip^^^ generated through almost de Sitter inflation 
is also 0{H). Using this relation, we have 



Thus we find that the CMB temperature fluctuation caused by the tensor perturbation 
directly probes the value of H during inflation, i.e. the energy scale of inflation. 
This can be discriminated from the scalar-type perturbation by looking at B-mode 
polarisation PUDIIII]. 

Another important observation is the non-Gaussianity of the temperature 
fluctuations fTZl [T3] . The non-Gaussianity is caused by the non-linear dynamics of 
cosmological perturbation. Once we have completely understood the evolution of density 
perturbation at late time, the remaining non-Gaussianity which is not accounted for 
should have their origin in the earlier universe. The results of WMAP 7 years indicate 
that the local-type non-Gaussianity parameter is given as f^L = 32 ± 21 at 68% 
confidence level [1]. For the Planck satellite, it is expected that the window of /nl 
is expected to be reduced to 0(10) |9]. 

Recently, non-Gaussianity of the primordial perturbation also has been studied by 



many authors [ZlElIiailSliailSllISllIZlIIHlIiniEOlEIlE^ 

EniElESlIMlEaESllMlE^EHlEBIinillllll^ The main reason 



for non-Gaussianity to attract much attention is the expectation to future observations 
mentioned above. These observations may bring us valuable information about the 
dynamics of inflation. Here in this short paper, we would like to clarify the difficulties 
in generating large non-Gaussianity from the non-linear dynamics during inflation. 

2. Basic feature of generation of non-Gaussianity of local-type. 
2.1. Local non-Gaussianity 

Primordial non-Gaussianity gives rise as non-trivial higher order correlation functions 
in primordial perturbations. In principle, various types of non-Gaussianity are possibly 
generated [251 EEl [27], but here we focus on the so-called local-type non-Gaussianity. 
Local-type non-Gaussianity is characterised by the existence of one-to-one local map 
between the physical curvature perturbation ({x) and the variable which follows the 
Gaussian statistics Cg{x) at respective spatial points. Namely, we can expand the 
curvature perturbation as 




(3) 



(4) 



and the Gaussian variable satisfies ordinary Gaussian statistics: 



(5) 
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Figure 1. Schematic picture to explain 5N formalism. 



In this case three point function can be characterised by a single non-hnear 
parameter /^vl [12] as 

(CkiCk.Cka) = ^^'Hki + k2 + k3)i?c(A;i, h), (6) 

with 

B^ih, h, h) = Pdki)Pdk2) + Pdk2)Pdh) + Pdh)Pdki) .(7) 

In general, other types of non-Gaussianity can be generated. However, when we focus 
on super-horizon dynamics, only local-type non-Gaussianity can be produced as is 
explained below. 



2.2. Super-horizon dynamics 

Here we just present an intuitive explanation about how the density perturbations evolve 
when the length scale is much larger than the Hubble scale. Super-horizon dynamics 
is locally described by the Friedmann-Robertson- Walker universe. We consider the 
evolution of the universe starting with an initial flat hypersurface at t = t*, on which the 
initial values of the infiaton field have a certain distribution. We evolve the spacetime 
until reaching the final surface at t = tp that is characterised by a specified energy 
density as shown in Fig. 1. Hence, the final surface is a uniform energy density surface 
by definition. As each horizon patch is causally disconnected from the others in the 
inflating universe, its evolution is determined as a local process. If initial conditions are 
completely given in each horizon patch, we can solve the evolution of its future. Here, 
for simplicity, we assume that the evolution of the averaged values of fields in each 
horizon patch is determined by the averaged values of initial data. This assumption 
is not necessarily true in general, but we have a good reason to assume so in most 
cases. Initial conditions for the smaller scale perturbations can affect the evolution of 
the averaged values of variables. However, as there are so many small scale degrees 
of freedom, the average of their effects will not be largely fluctuated except for rare 
situations. Of course, even if we can neglect the effect of fluctuations on small scales, 
this does not mean that the backreaction due to small scale degrees of freedom does not 
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modify the dynamics of the averaged values of variables. This backreaction effect arises 
at the second order perturbation or higher. Therefore the difference in the backreaction 
effect between different horizon patches is even higher order in perturbation. Thus, 
we usually neglect this effect for simplicity. Of course, in the situation hke preheating 
the sub-horizon scale perturbations evolve more rapidly than the averaged values, we 
definitely need to take the backreaction into account. However, even in that case, it will 
not be necessary to know all the details of small scale perturbations to understand the 
evolution of the averaged variables. 

Anyway, we neglect these effects originating from small scale perturbations here. 
Then, the e-folding number between the initial surface and the final one is completely 
determined by the initial averaged values of the variables. In the scalar field dominant 
universe, such initial conditions are specified by the values of the field components and 
their time derivatives, (j)\t^,x) and (p\t^,,x). We denote together with (j)\t^,x) 

by (j)\t^,,x). Then, the e-folding number between the initial surface and the final one, 

N{tF,U,(l)') = jj Hdt (8) 

is given as a function of (j)\t^:,x). Then, roughly speaking, the spatial metric on the 
final uniform energy density surface will be given by 

The curvature perturbation is in fact the perturbation in the above exponent 
N(tF,t^,(j)\t^:,x)), and hence we have 

atF)^6N{tF,<f)'{U,x)). (10) 

This formula is the heart of the 6N formalism [31 HI |5l |6l [71 [8] . 

In the slow roll case, the evolution equation for the scalar field can be approximated 
by a first order differential equation in t. Therefore the phase space of the initial 
conditions (f)\t^:,x) is reduced to a;). In this case, if there is only a single 

component, the trajectory in the field configuration space, which is one dimensional, 
is necessarily unique. Then, 5N{tp,(j)\t^,x)) does not depend on the choice of the 
energy density of the final surface. It can be set to a representative background value 
at the initial time. Namely, 

5N{tF, (p\U,x)) = 6N{U, (f)\u, x)). (11) 

Therefore nothing non-trivial happens during the super-horizon evolution for slow-roll 
single field inflation. Even if the slow-roll conditions are violated, single field inflation 
does not produce sizable SN for modes far beyond the horizon scale. This is because the 
trajectories in the configuration space continues to converge as the decaying perturbation 
mode gets irrelevant. In order to generate non-Gaussianity through the evolution on 
super horizon scales, it is therefore essential to consider multi-component infiaton field. 

Using the formula (10), we can express the curvature perturbations on the final 
surface in terms of the field perturbations on the initial surface as 

C(t J = 5N = N*5cPi + ]^NM5cPi + ■ ■ ■ , (12) 
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where the subscript * indicates the quantities evaluated at (j)^ = 4>^{t^:), and the subscript 
/ associated with N means the differentiation with respect to 



(13) 



where 4>\t) is the background trajectory. With the aid of Eq. (12), we can write the 
three point correlation function of ( a.t t = tp a.s 



(C(X1)C(X2)C(X3)) = N*N*jN*j, ((50i(xi)50,^(x2)50f (X3)> 

+ N;N*N*j,^ [(5</)^(xi)50,^(x2)<50f (x3)(50^(x3)> 



(3perms) 



+ 



(14) 



The first term contains only three 50s. If the initial fluctuations at t = t=K are Gaussian, 
this term vanishes. Depending on the model of inflation, early generation of non- 
Gaussianity at around horizon crossing is possible [22]. However, the observed almost 
scale invariant and slightly red spectrum of initial curvature perturbation strongly 
suggests that the inflation is likely to have been in the slow roll regime at around 
the horizon crossing. For the slow roll inflation, it is shown that the early production 
of non-Gaussianity is strongly suppressed [HI [151 [IS]- Therefore, it is well motivated 
to consider non-Gaussianity contained in the second term (or even higher order terms), 
generated by the non-linear evolution after the horizon crossing. 

Generation of non-Gaussianity can be classified into two classes. One is generation 
through the super horizon evolution during inflation. The other is generation at the end 
of or after inflation [^HlEniEIlElESlESllSZlllHlIll]. In either case the non-Gaussianity 
is of local type. Three point correlation function is characterised by one parameter f^L 
alone, and it is related to the e-folding number N{tF, 0^) as 



6 



NlN^N 



NL 



IJ 



(iVfiV^)2' 
neglecting higher order terms. 



(15) 



3. Non-Gaussianity produced at the end of or after inflation 

There are various models of inflation that produces large non-Gaussianity of local type 
at the end of or after inflation. Here we briefly mention curvaton type model [511 [52] as a 
typical example to explain the mechanism to generate large non-Gaussianity. Curvaton 
is another scalar field that is introduced to explain the origin of primordial curvature 
perturbation. During inflation, curvaton field is almost massless but it starts to behave 
as a massive field later when the Hubble rate decays down to the mass scale of the 
curvaton field. Then, the curvaton starts to oscillate around the minimum of the 
potential and eventually decays into radiation. At this point, the fluctuations in the 
curvaton a are transferred into curvature perturbations. For simplicity, we consider the 
case that the density perturbations are dominated by the fluctuations originating from 
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the curvaton. Neglecting a contribution of second order of Sp^ to (, which turns out to 
yield only /tvl = 0(1), we have [30] 

where r = Pa/ptot is the fraction of the energy density that the decay products of a 
takes. Here, we used the estimate pa ~ m^(cr + 6a)'^/2 with m being the mass of the 
curvaton. 

The amplitude of the power spectrum is observationally fixed as 

P^ = {[r6a/af) = 0(10"^°). (17) 

Therefore, if 6cr/a is not small, r can be as small as 10^^. On the other hand, the 
non-linear parameter is given by 

fNL = 0{{C)/{Cr)=0{l/r), (18) 

which can be as large as 10^. In this way, the production of large non-Gaussianity at 
the end of inflation is rather easily achieved. 

As a mechanism for getting large non-Gaussianity, one can also consider generation 
of perturbations during preheating phase [19], [50], SU HU HSj HH] HTJ.In this case, the 
background evolution is strongly coupled to the evolution of small scale fluctuations 
because the e-folding number changes depending on how fast the energy of coherent 
oscillation of fields decays into the small scale perturbations. While in the other models 
the transmutation of homogeneous isocurvature perturbations on each Hubble patch 
into adiabatic perturbations is determined just by solving the background evolution of 
various homogeneous and isotropic universe models. In this sense the mechanism to 
generate density perturbation during preheating stage is quite different from the other 
models. 



4. Non-Gaussianity produced during inflation 



Here we just present the formula for f^L in slow-roll inflation with canonical kinetic 
term. A systematic derivation for more general formula, including higher order 
correlators, is given in Ref. [21]. Assuming that the initial perturbations of fields are 
Gaussian in the following form 

{S<I>IM.) = + ^2)Psdh), (19) 

the formula for f^L is written in terms of the potential of the scalar field V as 
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NL 



Nf 



dN' Nj{N')Qi{N')Q\N')Q^{N'). (20) 

at = Np. To 



Here, we neglected non-Gaussianity in the relation between SN and 
define Ni and G*, it is convenient to introduce the propagator 



a; (AT, N') 



Texp / P{N")dN" 



N' 



(21) 
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where 



PUN) ^ 



— - + 

V 



y2 



(22) 



<t>=<l>{N) 

is the potential term that appears in the hnear perturbation equation when we use 
the e-folding number as a time coordinate. The indices associated with V represent 



differentiation with respect to scalar field. T in Eq. (21 ) means that the matrices are 
ordered in time when the exponential is expanded in power of Pj. Using this propagator, 
Ni and 9' are defined as 

N,{N)=N^{Nf)AI{Nf,N), 

e\N) = mN,N,)N^. (23) 
Ni{N) satisfies the equation of motion adjoint to the linear perturbation equation for 



_d_ 
dN 



N,{N) = -P^N,{N), 



with the boundary condition Nj{Np) 
same equation as (50* 
d 



(9iV(iV^,0O/50^')|<^=^(^, 



dN 



(24) 

e*(A^) satisfies the 
(25) 



and the boundary condition is given by G*(A^=i,) = N^{N^), where the index is raised by 
the inverse of the field space metric, which is assumed to be 5*-^ here. The three point 
interaction Q*^ is given by 



V 



+ 



(26) 



0=0(7V) 

In a naive sense, slow-roll conditions require that the potential of the infiaton is a smooth 
function of 0*. Hence higher order differentiations with respect to 0* are suppressed. 
More precisely, introducing a slow-roll parameter 

ivy' 



e = 



2 



(27) 



where iripi is set to unity, Pj and (5}^ are supposed to be of 0(e) and 0{e ' ), 
respectively. To obtain an almost scale-invariant spectrum, assuming this kind of scaling 
is natural although it is not strictly required [53]. In general, there might be a very 
massive isocurvature degree of freedom, which does not satisfy the above simple scaling. 
However, in this simple example, fluctuations in such a massive degree of freedom decay 
rapidly, and hence they can be safely neglected. 

Duration of the inflation is roughly estimated as 

f V \ ^ fHV 



N = 



(28) 



Hence, the exponent in the propagator A* is estimated to be 0(1). As it is the exponent, 
it is crucial whether the amplitude of this factor is slightly bigger or smaller than unity. 
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2 4 6 8 10 12 



Figure 2. A background trajectory in double inflation model. Taken from Ref. [2T] . 



In principle, therefore A*- itself can be much larger or much smaller than unity. The 
magnitudes of and 9-^ depend on the behaviour of the propagator, but the value of 
Ni ai N = Np is given by definition so as to satisfy 

V{^{Nf + 6N)) = V{^{Np)). (29) 
Expanding this equation up to first order, we obtain 

N,{Nf) = O {VV,/V'Vj) = 0{e-^/^). (30) 

If A* stays of 0(1), one can say that A^j and 6-^ also stay 0(e^^/^). Substituting these 



estimates into Eq. (20), we find that f^L = 0(e), as the general argument as given in 
Refs. [H Ha [IS] tells. 

However, if amplitude of A* largely deviates from unity, Ni and hence B* are also 
largely enhanced or suppressed, resulting in large enhancement or suppression of the 



integrand in Eq. (20). Then, the order of magnitude of /tvl is not necessarily suppressed 
like 0(e). Furthermore, slow roll conditions are violated, we have more chance to 
generate large non-Gaussianity. However, it is not so easy to construct a viable model 
which generates large non-Gaussianity from the super-horizon dynamics during infiation 
as we shall see below. 

5. Is there any model which produces large non-Gaussianity during 
infiation? 

5.1. Temporal slow-roll violation is not sufficient 

First we consider double infiation model ^1], which is a simple two-field model. The 
potential is given by 

V{<p,x) = \rnl<j)'' + \rnlx^. (31) 
There is no direct coupling between two fields, (j) and x, except for gravitational one. 
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Here we follow the analysis of Ref. [21]. We show the background trajectory for 
m^/m^ = 20 in Fig. 2 taken from the same reference. Since the mass of x-held is much 
larger, the field moves rapidly in x-direction first. Arriving at around the minimum in 
X-direction, the slow roll conditions are violated and the field oscillates several times. 
This violation of slow-roll conditions at the intermediate step is possible when the mass 
ratio is large enough. Otherwise, the trajectory becomes just a smooth curve and slow 
roll conditions are not violated in the middle. One may expect that this temporal 
violation of slow roll produces large non-Gaussianity. However, this is not the case. 

During the period when x-field is moving rapidly, 0-field stays almost constant 
because the mass of 0-field is small compared with the Hubble rate which is dominantly 
determined by the term m'^x^ the potential. Hence, the total e-folding number is 
simply given by the sum of e-foldings for two stages of inflation: 

N = Ar(<^) + N^^^ + xl (32) 



Applying the formula (15), we find 

Hence, as in the usual slow-roll case, the non-Gaussianity is suppressed by the slow-roll 
parameter e = 0(1/A^). 



5.2. Simple hybrid inflation does not give large non-Gaussianity 

As a second example, let's consider hybrid inflation model [HS] whose potential is given 
by 

^(0, X) = ^{X'- v^f + + ^0^^. (34) 

Usually (j) is the inflaton field and x stays at x = during inflation. The effective mass 
of the x-field is a function of given by 

ml = -Xv"" + g^(t)\ (35) 

assuming x = 0. The mass of x-held changes its sign at some point, which we denote 
by = Nc- Then, tachyonic instability in x-direction occurs, leading to the end of 
inflation. 

To analyse the dynamics of x-field at this critical point analytically, we approximate 
the mass of x-field as 

ml = -iJi\N -N,), (36) 

where /x^ is a parameter that depends on the choice of the model potential. 

If > iJ^, x-field is massive during (/)- field inflation. Thus, the fluctuation of 
X-field decays on large scales. If ji^ < if^, x field can stay nearly massless during 0-field 
inflation. In this case, the fluctuation of x-field may play an important role. As the 
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mass scale of x is small, we can apply the slow roll approximation to x-field, too. Then, 
we have 

3H^dNX = f^\N-N,)x. (37) 
This equation can be solved to obtain 

In ^ = [(AT - - (iV, - N,f] . (38) 

X* O-n 

Applying the formula (6/5) f^L ~ -^xx/^x' "^hich is valid when the fluctuation of x-field 
dominates, we have 

lfNL^^,{N,-N,). (39) 

Since we are restricted to the case with /i^ < H'^ and Np — cannot be very large, fjyi 
cannot be larger than 20 or so. If the fluctuations of 0-field also contribute, which is 
usually required in order to produce almost scale invariant spectrum, the value of fNL 
becomes even smaller. 

Furthermore, a more stringent constraint can be set by considering the following 
condition. If xi^c) ^ H is not satisfied, fluctuations at the horizon scale modes at 
that time determine on which side the field rolls down. Then, the tachyonic instability 
does not keep the coherence of x-held on the super horizon scales, leading to the so- 
called spinodal decomposition. As a result, long wavelength fluctuation is, roughly 
speaking, completely erased. Hence, there is no chance to generate large non-Gaussianity 
originating from the fluctuation of x-held. The condition xi^c) ^ H becomes 

^(Ar^-Ar.)2^1ng<ln|, (40) 

where we used the estimate xf ~ v. Using this inequality, the above estimate for f^i 
(39) is bounded by 



6. II 2 rupi 

5^"^~i7V3^"^' ^''^ 
where we used = 0{\v'^ /m^i) and |m^|0=o = Oi^Xv"^) > 0{H). The latter condition 

is required to avoid topological inflation happens at = x = 0. The factor ^| In 
is ~ 7 even if we reduce the energy scale of inflation ^JmplH to TeV scale. Here, as 
^ is also required, we conclude that f^L cannot be much larger than unity. 



5.3. Modular inflation 

In the preceding subsection we have observed that large mass scale for the isocurvature 
modes is necessary to generate large non-Gaussianity. However, the mass squared should 
not be large positive at the early stage of inflation when the interesting scales exit the 
horizon. This requirement is not easily compromised unless we invent a very artiflcial 
form of potential just to produce large non-Gaussianity. 
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Here we given one rather natural model which satisfies our requirement. This is 
the model that was proposed by Kadota and Stewart [56] ■ The basic model potential 
takes the form 

V = Vo- ml\^\^ + ^m|($3 + h.c.) + m||<l>|^ (42) 

where again we set rupi = 1. $ = (0 + ix)/A/2 is a complex scalar field. This potential is 
modified near $ = 0, the point with enhanced symmetry, due to loop correction of the 
Coleman- Weinberg type: Vioop = — In |$/<l>*|. This correction produces ring- 
shaped maximum around $ = 0, where topological infiation occurs. From the edges of 
the region of topological infiation, the field slowly rolls down the hill. 

Contour plot of the potential is given in Fig. 3. We focus on the trajectories which 
pass near the saddle point as shown in the same figure with a red curve. In this model 
mass squared in the x-direction is negative from the beginning. However, since it is the 
direction of approximate U{1) symmetry, the trajectories in any direction are almost 
equally preferred at this stage. As increases, the fiuctuation in x-direction is enhanced 
by the ratio of compared with its initial value. After the direction of trajectory 
changes, this fiuctuation in x-direction is transferred to the curvature perturbation. In 
this way the fiuctuations in x-direction can rather easily dominate the final density 
perturbations. Then, the amplitude of fiuctuations in x-direction stays almost constant 
near the ring-shaped top of potential, and therefore scale invariant spectrum is realised. 

Near the saddle point, the equation in x-direction can be approximated by 

3H^dNX = -^Ix, (43) 
with being the mass squared of X" field evaluated at the saddle point. Using the 
approximation =constant, this equation can be solved to give 

N=^\n^, (44) 
-m^ X* 
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where x* is the initial value of x- Applying the formula {Q/5)fNL ~ ^xxl-^x' obtain 
6 ml 

This can be large if ^ H^. However, to achieve this by natural form of potential, the 
potential minima as well as the saddles points should be located close to the origin in 
Planck unit. Then, inflation does not occur except for the top of the potential or near 
the saddle points. To sustain a sufficiently large e-folding number, inflation in these 
specific regions should last long. However, this is not allowed. If infiation lasts long 
near the top of the potential, which is ring-shaped, the motion in the radial direction 
becomes very slow. As a result, the amplitude of curvature perturbation originating 
from the fiuctuation in (j) direction dominates. For infiation to last long near the saddle 
point, the trajectory should pass very close to the saddle point. However, in such cases 
significant fraction of the universe will fall into the other side of the saddle point. As a 
result domain walls are formed, which leads to the problematic domain wall dominated 
universe. 

Even if we can avoid the domain wall formation (by considering models with higher 
symmetry), it is not easy to construct models which explain the naturalness of such a 
special trajectory. The probability distribution of the universes with different values of 
X* will be affected by the volume expansion factor during inflation, which is proportional 
to e^^. If 3i/^/|m^| is small as is requested for large non-Gaussianity, x* must be 
extremely small to gain a large e-folding there. Then, the probability of having such 
a small value of x* is extremely small. As long as we consider models which realize 
sufficiently large e-foldings for the natural choice of the fiducial value of x*i Inl will 
not largely exceed unity. Further variation of this model is possible by introducing loop 
correction at around the saddle points, but the final conclusion does not change much. 
We will report detailed analysis about this model in our future publication. 



6. Conclusion 



Multi-field infiation generates entropy perturbation as well as the adiabatic one. This 
entropy perturbation in general affects the evolution of the curvature perturbation even 
after the horizon crossing. Possible mechanisms of production of non-Gaussianity in 
the super-horizon regime can be classified into two cases. One is the production of 
non-Gaussianity during infiation and the other is that at the end of or after infiation. 
Although we could not give a general proof, it seems very difficult to produce large 
non-Gaussianity during infiation as far as we consider potential without fine-tuning, 
although there are several claims that suggest it possible [|j 

Here are several subtle issues. In this paper we claimed it difficult to construct a 
natural model in which large non-Gaussianity originating from the non-linear evolution 



I In the examples given in Refs. [40l |4T], it is not clear if we should say that the non-Gaussianity 
is generated during super-horizon evolution during inflation because distribution in field space stays 
Gaussian. In fact, the origin of non-Gaussianity is concentrated in the term neglected in Eq. (20). 
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is produced before inflation terminates. By contrast, generation of large non-Gaussianity 
at the end of or after inflation is rather easy. However, it may not be so trivial in general 
to identify when non-Guassianity is generated. We gave an example of double inflation 



in § Although we did not mention it in the text, in this model the non-linear 
parameter /tv^ as a function of tp becomes very large temporally when the background 
trajectory changes its direction, but this large fNL does not persist long. This means 
that the time dependence of /jvl is not always appropriate to read when non-Gaussianity 
is mainly generated. 

Another important point to note is the limitation of neglecting small scale 
perturbations. The small scale perturbations often become important for the field 
components that become massive during inflation. A typical example is the hybrid 



inflation model discussed in § 5.2 To take into account the effect of small 
scale perturbations, SN formalism based on the background evolution for spatially 
homogeneous spacetime is not sufficient. We need to extend 6N formalism incorporating 
collective variables that characterise statistical state such as the averaged magnitude of 
small scale fluctuation. 

Finally we should mention the naturalness of the initial conditions. It might be 
possible to obtain large non-Gaussianity by tuning the initial conditions, but it will not 
be realised in reality if they are extremely fine-tuned. As an example, we discussed 



Kadota-Stewart model in § |5.3[ However, once we start to discuss the likeliness of the 
chosen initial conditions, the long-standing measure problem arises. At the moment we 
are very far from a conclusive answer on this issue. 
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